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We derive a semiclassical trace formula for a symmetry reduced part of the spectrum in axially 
symmetric systems. The classical orbits that contribute are closed in (p, z,p p ,p z ) and havep^, = mh 
where m is the azimuthal quantum number. For these orbits vary with energy and almost 

^ ' never lie on periodic trajectories in the full phase space in contrast to the case of discrete symmetries. 

The transition from m = to m > is however not dramatic as our numerical results indicate, 
suggesting that contributing orbits occur in topologically equivalent families within which p v varies 
smoothly. 
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I. INTRODUCTION 



> 

Modern semiclassical theories deal with the duality of the quantum energy eigenspectrum and the classical spectrum 
of periodic orbit lengths and stabilities (![]. For integrable systems, the Poisson summation formula provides this 
connection and is equivalent to the EBK quantization scheme while for generic quantum systems, trace formulas 
deal with this duality and provide the only connection between quantum and classical mechanics. As an example, 
for hyperbolic systems the Gutzwiller trace formula expresses the density of quantum energy eigenvalues, d(E) = 
as: 

d{E) = d av (E) + -VV p cos (rS p - 7ra p r/2) (1.1) 

i nn pfri J\det(M;-I)\ 

where the summation over p extends over all primitive periodic orbits, T p is the corresponding time period, S p 
the action, a p the Maslov index associated with its invariant manifolds and M p the stability matrix arising from a 
linearization of the transverse flow. For systems where periodic orbits occur in families, appropriate modifications 
need to be made ||^] , though these may be system specific || . 

Quite often, the systems we encounter in nuclear, atomic or molecular physics have symmetries and one might be 
interested in the spectrum of a particular symmetry class. For example, billiard systems with reflection symmetry 
about the x and y axis have four classes of wavefunctions with a choice of Dirichlet or Neumann boundary condition 
on the symmetry lines. In these as well as in other systems with discrete symmetries, it is possible to construct the 
symmetry reduced Green's function and evaluate its trace [||. While the orbits that contribute are periodic in the 
reduced phase-space, they are not necessarily periodic in the full phase space. The end-points however are related by 
symmetry. It turns out nevertheless that when allowed to evolve in time, these orbits are eventually periodic in the 
full phase space. Thus trajectories contributing to any given symmetry class of the spectrum are either full periodic 
orbits or parts of these whenever the symmetries present are discrete in nature . 

The presence of continuous symmetries implies the existence of constants of motion. For example in systems with 
axial symmetry, the z-component of the angular momentum is conserved and quantum mechanically one may be 
interested in the spectrum of a particular m-subspace. A way to deal with this problem semiclassically is to work 
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with an effective two-dimensional potential V(p, z) + ft 2 m 2 / {2M p 2 ) where M is the mass of the particle and m is the 
azimuthal quantum number. Here V(p, z) is the actual potential in which the dynamics is executed and the additional 
term h 2 m 2 /(2Mp 2 ) arises from the conservation of l z [pU|]. The semi-classical trace formula then takes the form in 



Eq. (1.1) provided periodic orbits are isolated and unstable. 

An alternative procedure is to work with the dynamics in the full phase-space and it is desirable to understand the 
nature of the trace formula and the kind of classical orbits that contribute. We shall devote ourselves to this question 
through the rest of this paper and our results are in certain sense very different from those for discrete symmetry. 
The orbits that contribute necessarily have an angular momentum l x = mh and must close in p,p p and z,p z . Thus 
they are not necessarily periodic since they need not close in tp. Unlike discrete symmetries however, orbits that are 
not periodic are generally not parts of periodic orbits. We support these results with extensive numerical evidence. 

On completion of this work, we came across a paper by Creagh [|l0| that deals with this problem for arbitrary 
Abelian symmetries and rotational symmetry. Our approaches are different however. Creagh [ flO| takes recourse to 
group-theoretical methods in order to obtain the symmetry reduced Green's function and then uses a semi-classical 
approximation that ultimately involves classical trajectories possessing quantized values of these additional constants 
of motion. Our approach is based on the fact that for axially symmetric systems, the eigenfunction, ip(p, (p, z) = 
(2tt)~ 1 ^ 2 e lm!fi (j>(p, z) and we use this in the full Green's function to obtain the density of a given m-subspace. The 
restriction of constants of motion to their quantized values arises from a stationary phase condition in our case as 
also the fact that the trajectory must close in p,p p and z,p z . Thus, even though the final results are identical, the 
derivation in this paper provides an alternate viewpoint that is simple and transparent to those not familiar with the 
intricacies of extended phase space. 

Apart from a re-derivation of the trace formula and its numerical verification, this work emphasises the fact that 
orbits contributing to a symmetry reduced part of the spectrum are in general not eventually periodic in the full phase 
space whenever the symmetry is continuous. This is an important departure from the case of discrete symmetry. We 
also show that contributing orbits occur in topologically equivalent families within which p v varies smoothly. For 
successive m therefore, distinct orbits from the same family contribute and this is evident from the shift of peaks in 
the power spectrum of the quantum density. 



Ill 



The organization of the paper is as follows. In section |n|, we provide a derivation of the trace formula. Section 
deals with the inverse problem that motivate our numerical results which we present in section IV. Our conclusions 
are summarized in section 

II. FORMALISM 

Let us consider an axially symmetric three dimensional system described by the coordinates {ip, 77) where ip is the 
azimuthal angle about the axis of symmetry and 77 = (f^ - ^) represents the other two components in an orthogonal 
coordinate system. 

The basic definition of the full Green's function of such a system can be written as , 

n rn ' 

where 

^H^) = e -^P^^HV) (2.2) 

V Z7T 

In the above, both ip and (p are taken to be normalised and m is the azimuthal quantum number. Thus : 

G(*V, <pY; E) = J- £ exp[*m^" - <p>)) £ ^V^V) (2 3) 

ZTT L — ' z — ' £j — hj n m 



Multiplying both sides by exp(— ipp") and integrating over tp", one obtains, 

G(cp"n", tp'r)'; E) exp[-i(i(p" - <p')]d<p" (2.4) 



E 



E — E n 







where {E n ^} is the subset of eigenvalues for which the azimuthal quantum number equals p. Further integration 
over ip' on both sides gives, 
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e ^y^r — = f r r ^ ^ - ^ w. ^ 

E - E n ^ lit J J Q 

Eq. (|2.5| ) can be viewed as the symmetry reduced Green's function corresponding to the azimuthal quantum number 
fi. A semiclassical expression for this can be obtained by using the semiclassical approximation to the full Green's 
function 

G(<A", <ffrf; E) ~ £ \/|det tp'rf; E))\ exp[^S(E) - % -vv\ (2.6) 

^ ' cl.tr. 

which has contributions from all classical trajectories (cl.tr.) at energy E connecting points (ip" 7 rj") and (</?', rj'). The 
action S(E) is 

<p"v" 



S(r 1 " 1 r 1 \p,E) = f p.dq 
J u>'n' 



tp'T)' 

v" 

P n dr) + p v ((p" -tp' + 2ttN) 

■n' 

= ~S{r 1 ",r,',p tp ,E)+p v (<p + 2-KN) (2.7) 

where p v is a constant of motion due to the axial symmetry of the system, <p = tp" — ip' , N is the winding number 
and 5* is the reduced action f , p v drj. Lastly, v is the Maslov index which is determined by the caustics encountered 
along the trajectory while the matrix D is given by B] : 

D (ip'w\ipw-,E)= ^';' P ^ 

d(ip',r)',E) 

o 2 s d 2 s \ 

Lp'dr]" dip'dE \ 

^.o „ „o „ 1 



d 2 S 


d 2 S 


d 2 S 


dip' dip" 


d<~p' dr\" 


dip'dE 


a 2 s 


d 2 S 


d 2 S 


Or}' dip" 


dri'dr)" 


dii'dE 


d 2 S 


2 S 


2 S 


dEdip" 


OEdri" 


dW 



In order to arrive at a semi-classical form for the symmetry reduced Green's function, we insert Eq. fl2.6| ) in Eq. ( |2.5| ) 
and use a new set of variables ip = ip" — ip' and Tp = (ip" + p')/2 to obtain 

> = = -7^— Tvll^ ^cxp-wK/detf exp--— 2.9 



cl.tr. 

where we have used 

27T fill pTT pI'K 

/ dip' dp" = dp dp. (2.10) 
Jo J--K Jo 

We shall evaluate the integral over tp using the method of stationary phase in the limit h — > 0. The stationarity 
condition, 

pIQ 

--^ = (2.11) 

restricts trajectories to quantized values of p v (— dS/dip). We shall denote the stationary point by tp* such that 
p v (ip*) = [ifi. The integration yields 

y ^V^V) 1 \- [ 2 * Vldet^.T,",^;^)! 



£7--E„, M " (2^)3/2 7 ^ ^7y| 



riS 1 ivTr. .„ „„ 

x exp[- - — ] (2.12) 
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where U = d 2 S/dip 2 = dp^/dp evaluated at the stationary point ip* and v' 



(5' with P' = 1 if U is negative and 



zero otherwise. We shall postpone the Tp integration but for now we merely note that the symmetry of the system 
allows this to be evaluated exactly. 

The semi-classical symmetry reduced Green's function is thus expressed as a sum of contributions from classical 
trajectories at an energy E, connecting points 77" and 77' and having an angular momentum p v — /ih. We can now 
simplify the ratio of the amplitude determinants and express |det(D)| / \U\ as 



|det(J>)| 
\U\ 



det(D) 



2 s 



where 



D 



d^' dip" 



d 2 S 



dEdri" 



det(D) 




(2.13) 



(2.14) 



Details of this reduction can be found in the appendix. 

It is useful at this point to introduce a local co-ordinate system (772, 773) in the reduced (77) plane where 772 changes 
along the trajectory while 773 increases in the transverse direction and is zero on the orbit. It follows from the definition 
that 



d 2 S 
dEdrj' 2 



1 



d 2 S 
drj' 2 dE 



1 



and 



2 S 

so that the determinant ratio can be expressed as 



det(D) = 



1 1 



-( 



= = 



d 2 S 



d 2 s 

dri' 2 dr}'/ 



77) 



1 1 



I? 



(2.15) 



(2.16) 



(2.17) 



With these simplifications, we are now ready to evaluate th e tra ce of the symmetry reduced Green's function. 
Setting 77" = rf = 77 and integrating over 77 on both sides of Eq. (2.12), we get, 



1 



E n 



v ' cl.tr. u 



E 



2tt 



, 1 /T-57 AS iv'n 
drj — V \R exp — 

77 2 ft 2 



(2.18) 



where is the trace of the symmetry reduced Green's function corresponding to the azimuthal quantum number /z, 
& = § P-ndrj and 



R 



d 2 S 
dr)' 3 dri3 



(2.19) 



tp=tp*,ri'=Ti" 



Clearly, the only orbits that contribute to the trace are the ones that have p v — fjtfi and close in 77 but not necessarily 
in ip. Further, the 77 integration can be performed by the method of stationary phase and the stationarity condition 
then picks only that subset of orbits for which p v i — p n » at 77' = 77" = 77. We shall refer to such trajectories in the 
full 31? dynamics as quasi-periodic trajectories (q.p.t.) since they will be periodic in the reduced dynamics of the 
77-motion. In what follows such periodic orbits in the reduced system will be referred to as reduced periodic orbits 
(r.p.o.). 

The action S in the neighbourhood of a r.p.o. can thus be written to the second order as 5(773, £7) = 5(7/3 = 



0;E) + \W where 



W 



d 2 S 



d 2 s 



d 2 s 



dV3 d V3 drjl dr]' 3 dr]' 3 d^ 3 



i 3 =13= 



(2.20) 
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On performing the 773 integration by the method of stationary phase, Eq. ( [2.18 ) reduces to 



9»{E) 



2m% 



E 

q.p.t 



dtp (f 

J m 



R 




W 





exp[-5(£) 



7T 

l 2 a] 



(2.21) 



where a = v 1 + /3 with /3 = 1 if W is negative and zero otherwise. As in the case without symmetry, a is the Maslov 
index of the stable or unstable manifold and is an invariant of the reduced periodic orbit [fy]j|] . 

While we do not explicitly show this, the factor involving the second derivatives of S is related to linearized dynamics 
in the neighbourhood of the r.p.o. and is expressed as M: 



\ 



3 2 S 



dr, 3 dn 3 



^3 =J ?3= 



d 2 S 

a I) a T, 



' dru dv' 3 



d*S 
drudri s 



V|det(M-/)| 



(2.22) 



where the 2x2 matrix M is the stability matrix describing the dynamics in the linearized neighbourhood of the r.p.o. 
and I is the unit matrix. Note that tr is independent of the position 7/2 along the periodic orbit and so are S(E) and 
y/\ det(M — I)\. The integration along 772 thus yields 



dm 

m 



(2.23) 



where To is the period of the primitive r.p.o., or, in terms of the time, T q , p , t . required to span the q.p.t. in the full 
3D -dynamics, To = T q . p .t./No where Nq is the number of primitive r.p.o. contained in the quasi periodic trajectory. 
It is instructive however to deal with the 772, integrations together and this yields the "volume" of initial conditions 
for which the various stationarity conditions are satisfied. Thus, 



2,T 



dtp<b^ = W 



(2.24) 



where Wo = $T = 27rTo = 2TrT q . p .t./NQ. Note that when the q.p.t. is also periodic and has an N\ fold symmetry, one 
can interpret $ = 2n/Ni and T = T q , p , t /r = T p when viewed in the full phase space. Here r is t he repetition number 
of the periodic orbit in the full phase space and rN\ = Nq. With these clarifications, Eq. (2.21) becomes, 



1 



2mfi 



E 

q.p.t 



2ttTo 



V|det(M- J)| 



exp [-<?(£) 



(2.25) 



This is the main result of this section. The sum involves quasi-periodic trajectories at an energy E and with p v 
and is identical to the result obtained by Creagh [n0| . 



It is instructive however to examine the significance of Eq. ( 2.11 ) in greater detail. In performing the tp integration 
for each orbit, it is implicit that the action changes smoothly as the orbit (characterized by tp and hence p v ) is varied. 
In other words, orbits occur in topologically equivalent families fl2]| within which the action varies smoothly and the 
stationary phase condition of Eq. ( [2.11 ) picks out one orbit from this family. By changing /i, another orbit from the 
same family contributes and this has a slightly different value of S. The range of fi over which a family contributes 
depends on its extent which in turn is decided by the potential V(p,z). These ideas are put on a more concrete 
footing for axially symmetric cavities in the following sections. 

It is easy to check that the trace of the full Green's function is recovered by summing over /i as follows : 



AT ** 



N 
1 1 

ih 2-7T 



TV 



E/^E 



27TTo 



^ ^|det(M-/)| 



(2.26) 
(2.27) 



where we have used the Poisson summation formula and substituted Eq. ( 2.25| ) for g M (E) . The integral when evaluated 
by the method of stationary phase picks up only those trajectories for which dS/dp v + 2nN = or tp* = 2irN. Thus 
the only orbits that contribute to the full spectrum are the ones which are periodic in the full phase space. The final 
expression is 
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9(E) 



(2tt/JVi)T p 



<9p ¥ 



9y> 



exp[-S p (£) 



1 ' 1 
2^ 



(2.28) 



with a' = a — 5 where S = 1 or depending on whether d 2 S/dp^ is positive or negative. The result is identical to 



ift(27rift)V2 ^ v /|det(M-7)| 

that obtained by Creagh & Littlejohn p.p^[. 

For the sake of completeness, it is necessary to mention that there is yet another contribution to the trace that we 
have so far neglected. It arises form the zero-length orbits and gives rise to the average density of states, d£ v (E) and 
to the leading order this is given by H,10| : 



1 



dpdq h 5(p v — ph)S(E — H(p, q)) 



(2.29) 



Eq. ( 2.25| ) gives the oscillatory part of the density, d% sc (E) = — ^ lim e _ +0 + 3 g p (E + ie) where 9 denotes the imaginary 
part. In the following sections we shall illustrate our results for axially symmetric cavities and look at contributions 
of q.p.t. in d% sc (E) as well as the p, dependence of d£ v (E). 

Though we have restricted ourselves to the case of axial symmetry, the formalism is very general and can be adapted 
to other cases of continuous symmetry. For example in case of rotational symmetry, one can express the wave function 
as : 



(2.30) 



where {Y, m (9, tp)} are the spherical harmonics. This can be used in Eq. (2.1) to obtain the trace of the Green's 
function for a given value of the quantum numbers I and m. 



III. AXIALLY SYMMETRIC CAVITY 



We now specialize to the case of a particle moving freely inside an axially symmetric cavity and under going specular 
reflections on collision with the wall. As in most other systems, it is easier to verify the duality in Eq. ( 2.25 ) starting 
with quantum levels and obtaining the length spectrum of periodic orbits. This is easily achieved by computing the 
power spectrum of the quantum density, d fl (v r E) : 



S"(i) 



dkd"(k)exp{ikl} 



J dk(d% v + dUeMikl} 



(3.1) 
(3.2) 



where \JE = k and eP(fc) = 2kd p (k 2 ) is the quantum density for a given value of /i. The average part of the quantum 
density gives rise to a peak at zero in S^Q) while the oscillatory part contributes peaks at the length of quasi-periodic 
orbits with p v — fih. 

What is perhaps not apparent is the fact that for [i ^ 0, the reduced periodic orbits that contribute to d^ sc vary 
with energy whenever they are not periodic (i.e. tp 0). These peaks are therefore expected to shift in the power 
spectrum with the window over which it is computed. 

Consider for example an orbit in the equatorial plane [z — 0) where the particle reflects specularly inside a circular 
billiard of radius, R. It is easy to see that two successive points on the caustic [[l5| (one reflection between them) 
satisfy the conditions for quasi-periodicity since p' = p" , z' = z" — 0, p' p — p" p and p' z = p" z = 0. The action on this 
q.p.t is thus : Sj = \f~Elj where lj = 2Rsvn(p/2) where p is the difference in angle, p" — p between these points on 
the caustic. By moving away from the caustic to a point along the trajectory, there exists another point separated by 
a collision where conditions for quasi-periodicity are satisfied (i.e. p' — p" and p' p = p" p ). It is also easy to check that 
along this quasi-periodic orbit, p is conserved while p changes by moving from the c aust ic in the radial direction to 
another orbit of the same topological family. The stationary phase condition in Eq. ( [2.11 ) is thus : 



ERcos{p/2) = ph 



(3.3) 



or equivalently : 



p* = 2 cos- 1 (ph/RVE) 



(3.4) 



G 



so that 




Sj = 2B.VE\ 1 - ^rr; - 2ph cos" 1 (uti/RVE) (3.5) 



in Eq. (2.25). For p, ^ 0, the orbit varies with energy continuously and the system is thus inhomogeneous. At a 
given (non-zero) value of p, the length of the contributing orbit increases with energy and thus should give rise to 
broadened peaks in the power spectrum, S' A1 (/). 

Note that the stationarity condition in Eq. ([O]) picks up one orbit for each p from this topologically equivalent 
family of 1-bounce orbits within which p v varies smoothly. In spite of the complicated dependence of S on E and p, 
it is evident that with increasing p, peaks in the power spectrum, 5 M (Z), shift to the left. Moreover the broadening of 
peaks is expected to increase with p due to larger spread in the length of contributing orbits. 

An interesting consequence of this analysis is the fact that at almost all energies this quasi-periodic trajectory 
is not eventually periodic in the full phase space whenever p ^ 0. This follows from the fact that periodic orbits 
in the circle billiard form a set of measure zero. This conclusion also holds in the more general case. For ^ ^ 0, 
whenever a reduced periodic orbit corresponds to a periodic orbit in the full phase space with p v = ph, the action, 
Sj = \J~E lj — 2TrNjfj,h where Nj is the winding number of the orbit. However, on changing the energy, E, this orbit 
no longer has the quantized value of p v and hence cannot contribute. Similar arguments hold for q.p.t's which form 
part of periodic orbits in the full phase space. Thus at any energy, most orbits that contribute to d% sc do not lie on 
a trajectory that is periodic in the full phase space. This phenomenon is in sharp contrast to the case of discrete 
symmetries where all quasi-periodic trajectories are eventually periodic in the full phase space. The only exception is 
the trivial case, p = when all reduced periodic orbits lie on periodic trajectories in the full phase space. 

Though we have restricted ourselves to a simple orbit in the preceding analysis, the conclusions are in fact more 
general as the numerical results of the following section indicate. 



We now return to an evaluation of the average density of states, d% v (E) for cavities where Eq. (2.29) becomes : 

If p 2 

d% v (E) = J dp p dp z dp v dzdpdip % 5{p ip - ph)5{E - p p — p z — -|) (3.6) 

for a particle of mass, M = 1/2. The integrals are straightforward to evaluate and the final result is expressed as : 

= -7-tz-A! (3.7) 
47m 

where 

A' = dz dp. (3.8) 

Note that , z^ ax , p^ in depend on \i and are dictated by the centrifugal barrier while p m ax depends on z through 
the shape of the cavity. For small \p\ j^fE however, A' ~ A — ZqRq where Zq is the extent of the cavity along the 
Z-axis and R — \p\ h/y/E. Thus 



where 



A = [ m "* dz ( P dp (3.10) 

Jz Q min h 

is the area of the pZ plane. Consequently, the average integrated density of states, Ng v (E) — d^ v {E')dE' is : 

where for convenience we have set % = 1. 

Note that the first term in the expression for N£ V (E) is identical to that for a 2— dimensional billiard of area A. 



We have not evaluated the standar d peri meter term here which is of the same order as the second term of Eq. (3.11) 



though as we shall show later, Eq. (3.11) predicts the correct p dependence for cavities. 
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IV. SOME NUMERICAL RESULTS 



We present here some numerical results on an axially symmetric cavity described by 

R(9) = ^{l + a n P n (cos 9}) (4.1) 

where Rq is the radius of a sphere having the same volume as the deformed cavity, A is a volume preserving factor, 
a n is the strength of the deformation and {P n } are the Legendre polynomials. R(6) is the distance from the origin 
to a point on the surface of the cavity, which due to the axial symmetry depends only on the angle, 9 measured from 
the Z-axis. For our calculations, we choose a P2 (n = 2) deformed cavity with Rq = 1, A = 1.008087 and = 0.2. 
For convenience, we choose the mass, M = 1/2 and % = 1. Details of the procedure used in obtaining the quantum 
energy eigenvalues can be found in [^6| . 

We have made no attempt towards finding periodic orbits systematically for this system both because there is no 
obvious symbolic dynamics that one can use and also since we are interested in the inverse problem. While orbits 
in the XY plane are well known, we have determined periodic orbits in YZ plane (p v = 0) using the orbit length 
extremization technique. We believe, the list we have compiled is complete for lengths, I < 15. Fig. (1) shows a plot 
of two periodic orbits which also satisfy the condition p v > — p n n at points separated by half their total lengths. 

In Fig. (2-4) we show plots of the power spectrum, S /i (Z) for p = (0, 1), (2,4), (6,8). The small arrows mark the 
position of full periodic orbits in the YZ plane and the orbit along the diameter in the XY plane, while the larger 
arrows mark those orbits that do not close in ip. All of these h&vep v = and hence orbits which do not close in tp form 
part of periodic orbits. The peaks in the power spectrum for p, = agree very well with the lengths of quasi-periodic 
orbits marked by arrows though there is no distinct peak at 4.7615 or 9.5230. These correspond to the periodic orbit 
along the Z-axis and its repetition. Unlike other orbits which occur in a 1-parameter family, this orbit is isolated and 
hence has a much smaller contribution. We shall return to this later while discussing states belonging to particular 
parity class. 

The /i = 1 power spectrum is also described well by these orbits which have p v = 0. The peaks shift to the left 
only slightly indicating that the topology of the orbit does not chang e wh ile going from p v = to a non-zero value of 



P(p. The phenomenon can be understood from the analysis in section III where we have explicitly described the orbit 
selection process for the simplest orbit in the XY plane. The only significant change occurs at 10.38 where the peak 
height decreases appreciably. 

This process continues for larger values of p, and can be seen from the power spectrum. The peaks gradually shift 
to the left of the arrow and nothing dramatic happens except that at p, = 4, there is no trace of the orbit at I — 10.38. 
There is yet another feature that is readily seen by comparing the power spectra of p = or 1 and p — 8. The peaks 
broaden significantly with increasing p even though there is hardly any change in the range of eigenvalues considered 
for evaluating the power spectrum. This is due to the fact that the length of the orbit (belonging to a topologically 
equivalent family) varies with energy and its spread around the mean value increases with the value of p. 

We emphasize this aspect in Fig. (5) where we plot the power spectrum for p = 26 for two different windows of 
energy. The one centred at the higher energy has peaks significantly to the right of the other, indicating that the 
lengths of contributing orbits increase with energy. 

We also compare the power spectrum of p = 1 and p = 26 in Fig. (6) and find that for the simplest (topological) 
orbit family (connecting successive points on the caustic) the peaks for p = 26 are broader and to the left of the 
corresponding peak for p, = 1. 

These findings corroborate our results of the previous sections and we now investigate the power spectrum of 
even/odd parity states for a particular value of p. Contribution to the semiclassical symmetry reduced even/odd 
Green's function come from two sources. Apart from the usual orbits with p v = ph connecting (p',z',(p') 
and (p",z",tp") there are additional contributions from orbits having p v = ph connecting points (p',z',ip') and 
(p" , —z", <p" + 7r) weighted appropriately by the character of the symmetry group |5j. Alternately, one can exploit the 
additional tp symmetry and express the symmetry reduced even/odd (denoted respectively as +/-) Green's function 
as : 



E 



^(p,z)^ n f(p,z) = 1 
E - Eft u p 4^ 



r 2n r 2n 



2tt i-2-k 

G(tp ,p, z; <p\ p, z; E) exp[— ip((p — ip)]dip dtp 



Jo 



±(-1)" / / G(Lp",p,-z;<p',p,z-,E)eM-w(<P"-<p')W'd<p' 
Jo Jo 



(4.2) 



The trace of its semiclassical version thus involves quasi-periodic trajectories with p v = \xh and trajectories connecting 
(p, z) with (p, —z) having p v = pjh, p{ — —p\ where / and i refer respectively to the initial and final points. 
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Thus we expect to see both full and half q.p.t.'s in the power spectrum of even parity and odd parity sequences. In 
particular, the half orbit along the Z-axis, will have a lower damping (due to its instability, the damping is exponential 
with the length of the orbit) and hence could be visible in the power spectrum. This is indeed obvious in Fig. (7) 
where a peak appears at I — 2.38. In fact, a closer observation also reveals a peak at I — 4.76. 

Significantly, orbits belonging to the XY plane do not have any other symmetry related orbit and hence have no 
additional peak at half the length. 

Finally, we study the [i dependence of the average density by plotting kffl — kffl as a function of the level number, 



N in Fig. (8). Here fcj^- = \/E%j- An approximate expression for — k*P can be obtained using Eq. (3.11) and is 



Tjjy y iJjy. .ni.ll apjJlUAllUdlC CJVJJl G^niUll 1U1 rujy fl/^y 

given by : 



+ ^ + (4. 3) 



Therefore for large N, fcjy — fcjy ~ \jA^o_ an( j f or ^ e cav ity that we have considered, this is 1.38 |/z|. 

Fig. (8) is a plot of k^p — kff for [i = 1 (bottom curve) to fi = 6 (top curve) . Each curve is approximately constant 
and the mean separation between curves is 1.38 \fi\ as predicted. 

V. CONCLUSIONS 

In the previous sections, we have studied the semiclassics of symmetry reduced spectra for the case of axial symmetry. 
There are two approaches to the problem. The first exploits the conservation of l z to reduce the system to 2-dimensions. 
Orbits then move around in an effective potential that includes the centrifugal term h 2 m 2 /(2M p 2 ). In comparison, 
the approach that we adopt here uses the dynamics in the full phase space and the trace formula then takes a different 
structure altogether. Our aim has been to understand the nature of the classical trajectories that contribute and also 
point out the differences with the case of discrete symmetry. 

Most of this paper deals with a derivation of the trace formula for the spectrum with a given value, p of the 



azimuthal quantum number, m. Though this has been studied by Creagh 10 earlier, we provide an independent 
derivation here using only the structure of the wavefunction to achieve symmetry reduction. Finally, though we have 
worked with axially symmetric systems, the method can in principle be extended to a large class of systems with 
continuous symmetry. 

Our main results are the following : 

• The average density, d£ v (E) = £ - lj ^r^ for small \p\ /y/E so that k$ - k$ ~ ^ where N denotes the 

level number and k^ = \/En- 

• The only classical trajectories that contribute to the oscillatory part of the quantum spectrum, d% sc (E) are those 
which close in (p, z,p p ,p z ) and for which p v — ph. Thus they need not be closed in the full phase space though they 
may be eventually periodic when allowed to evolve. These orbits are however periodic in the reduced phase space 

{P,Z,Pp,Pz)- 

• For p, = 0, all reduced periodic orbits are identical to or part of periodic orbits in the full phase space. 

• For fi 7^ 0, almost all orbits that contribute to d^(E) at any given energy are non-periodic in the full phase space 
even when they are allowed to evolve in time. Peaks in the power spectrum thus shift with the window over which it 
is computed and the width of peaks increases with p. 

Our numerical results support these observations and allow us to state the following : 

• for smooth potentials (billiards with smooth boundaries) , topologically equivalent orbits contribute at successive 
values of p as can be seen from the gradual shift in the position of peaks in the power spectrum. This is implicit in 
the selection criterion in Eq. ( [2.1l| ) since each of the orbits selected necessarily occurs in a (topological) family within 
which p v varies smoothly. 

The transition from fj, — to fi ^ is therefore smooth and this might have some significance in quantizing the 
spectrum for fj, = 1 using information about contributing orbits at /i = 0. 
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APPENDIX A: 



We outline here the steps involved in the reduction of the amplitude determinants leading to Eq. ( 2.13] ). 
The equation for the reduced action 

S(r,",r,',p v , E) = S(r)",r)', <p, E) - p v p (Al) 

is a Legendre transformation from the variables (rj" , rj , ip, E) to the variables (77", rj' ,p v , E) and leads to the relations 

dS__dS_ dS__dS_ dS_ _ dS_ dS_ _ dS_ 

; Wl ~ drf!' dE~dE ; lfy~ Pv ' dp~^~ 

The following expressions involving the second derivatives of the action can now be derived : 



(A2) 



d 2 s 


d 2 S 


d 2 S d 2 S d 2 S 


drfidrj'J 




dppdrii dp^drj" dp 2 


d 2 S 


d 2 S 


d 2 S d 2 S d 2 S 


dr&dE 


d^dE 


dp v drl i dp v dE dp 2 


d 2 S 


d 2 S 


d 2 S d 2 S d 2 S 


drfidE 


di)'!dE 


dp v dr}" dp v dE dp 2 


d 2 S 


d 2 s 


d 2 S d 2 S d 2 S 


dE 2 


~ dE 2 + 


dp v dE dp v 8E dp 2 




d 2 s _ 

dip 2 


\dp 2 J 




d 2 s 


d 2 S d 2 S 




d?f i dp 


dp^dr/l dp 2 




d 2 S 


d 2 S d 2 S 




drfldip 


dpipdr]" dp 2 




d 2 S 
dtpdE 


d 2 S d 2 S 
dp v dE dp 2 



We shall derive Eq. ( |A3a| ) as an illustration. 

8 2 S 



{rf',rf,p v ,E) 



Also, 



d_ I dS 

d ( dS 

d 2 S d 2 S d Pv 
drj^drj'J dp v drj" dr( t 



d 2 s d fas „ , 



(A3a) 



(A3b) 



(A3c) 



(A3d) 



(A3e) 



(A3f) 



(A3g) 



(A3h) 
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so that 



Thus, 



d 2 S dip 
dipdrt dp v 



d 2 S d 2 S 
dptpdr)'- dtp 2 ' 



d 2 S 



d 2 S 



d 2 S d 2 S d 2 S 



dr^dr)'! drj^d-q 1 ' dp v di]" dp^dr/^ dip 2 



(A4) 



(A5) 



We a re n ow in a position to deal with the ratio of determinants. By interchanging rows and columns, the determinant 
of Eq. (|2.8|) can be expressed as 



det(D) = dct 



d 2 s 


a 2 s 


d 2 S 


dr\' dr\" 


dr)>OE 


drj'dip" 


d 2 S 


d 2 S 


a 2 s 


dEdri" 


dW 


dEdip" 


d 2 S 


a 2 s 


d 2 S 


dip'dri" 


dip'dE 


dip' dip" 



d 2 s 



dtp'dp 



-det (£>-£) 



(A6) 



where 



and 



E= 



D 



d 2 s \ 

drj'dip" 
d 2 S 
\ dEd<p" ) 

( 8 2 S ' 

ydip'dip", 



d z S 



d z S 
dEdri' 



d 2 S 
drfdE 
d 2 S 
dW 



(A7) 



/ d 2 s 

\dip'dp" 



( d 2 S 
y dp' drj" 



d 2 S 



d 2 S 



d 2 S 



irj'dip" dip'drj 1 
d 2 S d 2 S 



dp'dE 

d 2 S d 2 S 
drj'dp" dp'dE 
d 2 S d 2 S 



dEdp" dip'd-q" dEdip" df'dE 



Using relations (A3a)-(A3h), it follows that 

, . - , det D 
det(D) = -gg- 

dtp 1 dip" 



det(D -E)= det 



d A S 
dr]'drj" 

d 2 S 
dEdri" 




(A8) 
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FIG. 1. Two periodic orbits in the YZ plane on which points separated by half the total length satisfy the conditions, 
z' = z",p' — p" and p z i — p z ",p p ' = p p ". Their respective lengths are 11.532715 and 17.856224. 



FIG. 2. Power spectrum, 5 M (i) of p — & 1 states. The shorter set of arrows mark the lengths of full periodic orbits in the 
Y Z plane while the longer arrows are quasi-periodic trajectories in the YZ plane that do not close in ip. All trajectories have 
p v = 0. 



FIG. 3. Power spectrum, 5 M (/) of p — 2 & 4 states. The arrows mark the length of orbits with p v = as in Fig. (2). 



FIG. 4. Power spectrum, S^(l) of p — 6 & 8 states. The arrows mark the length of orbits with p v = as in Fig. (2). 



FIG. 5. Power spectrum of p — 26 states for two different ranges of energy. The window centred at higher energy has peaks 
shifted to the right. 



FIG. 6. A comparison of the power spectrum of p = 1 and 26 states. The p — 26 spectrum gives rise to broader peaks due 
to a larger spread in the lengths of contributing orbits. The window is approximately over the same range of energy in both 
cases. 



FIG. 7. Power spectrum of p = 0, even parity states. The smaller arrow marks the length of the half-orbit along the Z-axis. 



FIG. 8. A plot of k£ — &n as a function of the level number, N. The curves from bottom to top are (in order of increasing 
p) for p — 1 to p — 6. 
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FIG. 1. Two periodic orbits in the Y Z plane on which points separated by half the total length 
satisfy the conditions, z' = z",p' = p" and p z i = p z ",p p i = Pp»- Their respective lengths are 
11.532715 and 17.856224. 
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FIG. 2. Power spectrum, S^(Z) of \i = & 1 states. The shorter set of arrows mark the lengths 
of full periodic orbits in the YZ plane while the longer arrows are quasi-periodic trajectories in 
the YZ plane that do not close in tp. All trajectories have p v = 0. 
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FIG. 3. Power spectrum, 8^(1) of fj, = 2 & 4 states. The arrows mark the length of orbits with 
p v = as in Fig. (2). 
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FIG. 4. Power spectrum, S^ty) of = 6 & 8 states. The arrows mark the length of orbits with 
p v = as in Fig. (2). 
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FIG. 5. Power spectrum of \i = 26 states for two different ranges of energy. The window centred 
at higher energy has peaks shifted to the right. 
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FIG. 6. A comparison of the power spectrum of \i = 1 and 26 states. The \i = 26 spectrum 
gives rise to broader peaks due to a larger spread in the lengths of contributing orbits. The window 
is approximately over the same range of energy in both cases. 
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FIG. 7. Power spectrum of \i = 0, even parity states. The smaller arrow marks the length of 
the half-orbit along the Z-axis. 
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FIG. 8. A plot of k$ - k$ as a function of the level number, N. The curves from bottom to 



top are (in order of increasing jj) for \i = 1 to \i = 6. 
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